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By Sidney M. Harmon 


SUMMARY 

A method utilizing source singularities is presented for obtaining 
the linearized dowitwash field due to lifting vings of infinitesimal 
thickness at subsonic and supersonic speeds. The distribution function 
for the source singularities is specified by the loading on the ving. 
The method is applied to derive generalized formulas for the down-wash 
field due to uniformly loaded swept and rectangular -wings at subsonic 
and supersonic speeds. The utilization of these formulas to ob-tain the 
do-wnwash due to wings of arbitrary loading is indicated. An example of 
the procedure is given in which specific formulas are derived for the 
downwash field due to a rectangular -wing at supersonic speeds for a 
uniform loading and for a linear chordwise variation in loading. 


INTRODUCTION 

Several methods based on linearized theory are available to obtain 
the downwash field due to lifting surfaces at subsonic and supersonic 
speeds for use, as an example, in stability calcxlLations. The calcu¬ 
lation of the downwash field at subsonic speeds has relied almost 
exclusively on Prandtl's lifting-line theory, which is based on the 
concept of a horseshoe vortex (for example, reference 1). Present 
methods for calculating the do-wn-wash field at supersonic speeds are 
those utilizing conical flows (reference 2), potential doublets (refer¬ 
ence 5), vortices (references 4 and 5), and pressure doublets (refer¬ 
ences 6 and 7). The integrations required of the foregoing vortex 
or doublet singularities or in the conical-flow method in order to 
obtain exact solutions of the linearized equations for lifting surfaces 
have generally been found to be difficult} therefore practice has 
usually had recourse to approximate methods based on lifting-line 
theories (references 5 and 7). 

The present report, prepared at the MCA Le-wis laboratory, indi¬ 
cates a method that is intended to facili-bate the computations for 
obtaining the exact linearized downwash field due to lifting surfaces 
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at subsonic and supersonic speeds. The method utilizes source singu¬ 
larities with the distribution function specified by the loading on 
the surface. The method is applied herein to derive fonnulas for the 
downwash field due to uniformly loaded swept and rectangular wings of 
infinitesimal thickness at subsonic and supersonic speeds. The utili¬ 
zation of these formulas to obtain the downwash field due to wings of 
arbitrary loading, by means of the correspondence relations presented 
in reference 8, is then indicated. An exan^le of the procedure is 
given in which specific formulas are derived for the rectangular wing 
at supersonic speeds for a uniform loading and for a linear chordwise 
variation in loading. 


SYMBOLS 


The following symbols are used in this report: 

A,B,Ci,C 2 refer to regions bounded by foremost Mach aft cone 

^ corners of rectangular wing at supersonic speed 

2^*** 5 (plan view of regions in fig. 3) 

a constant used to describe prescribed linear chordwise 

variation in loading 

g := - 1 (also used to refer to Mach cone region as 

indicated in fig. 3) 

■R. = Vl - 


arbitrary constant 
chord for rectangular wing 


h wing semispan 

I refers to integral and solution expressed by equa- 

^ tions (25b) and (25e), respectively 

I refers to integral and solution expressed by equa- 

^ tions (25c) and (25g), respectively 


K 


constant in equation (l) (K = 
K = — at supersonic speeds) 

7t 


1 

2jt 


at subsonic speeds; 
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M free-streaiTi Mach niimber 


cot A y where A is angle of sweep of wing leading 
edge (fig. 1) 


r = 

)^-B^ [iy-Tl)“+z3 

r’ 


u,v,w 

disturbance velocities of fluid in flow field aG-ong 

X-, Y-) Sell'S, z-axes, respectively (fig. 1) 

^0 

value of u on upper surface of wing for uniform 
prescribed loading 

w 

refers to downwash due to left half of wing 

^0 

refers to downwash due to uniform prescribed loading 
on wing 


refers to downwash due to linear chordwise variation 
in prescribed loading on wing 


contribution to downwash of continuous portion of 
leading edge for semi-infinite oblique wing with 
uniform loading shown in figure 1; also represents 
downwash due to semi-infinite line source originating 
at leading edge of center section 


contribution to downwash of discontinuity of leading 
edge at origin in semi-infinite oblique wing with 
uniform loading shown in figure 1 


downwash due to term in series formed by expressing u 
on wing as function of x, y, and uq 


contribution to downwash of streamwise side edge for 
semi-infinite oblique wing with uniform loading 
shown in figure 1; also represents downwash due to 
semi-infinite streamwise line source originating 
at leading edge of center section 


represents downwash due to semi-infinite line source 
originating at (xQ,yQ) on right and left halves of 

wing, respectively 
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x,y,z 


x',y',z’ 


Xq^Yo 


rectangular coordinates with origin at leading edge 
of center section (fig. 1) 

oblique coordinates related to rectangular coordinates 
according to equations (5) 

rectangular coordinates indicating origin of semi- 
infinite line source 


) contribution to downwash of wing cut-off at trailing 

^ edge for uniform loading 

£ infinitesimal distance in y-direction across side 

edge 

auxiliary variables used to replace x, y, z, x% 
y*, and z', respectively 

= V|b^ 1 (y^+z^) 

e = VIb^i (y'2+z'2) 

i upper limit for integral in equation (7) (See dis- 

cussion following equation (7).) 

T region of integration (See discussion following 

equations (l) and (3) •) 

function representing solution of linearized partial 
differential equation (2) 


Subscripts: 


A,B,C ,C refer to corresponding regions indicated in figure 3 

1 2 


* * *^5 
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Cj. refers to chord at center section of swept wing 

c^ refers to chord at tip section of swept wing 

L,S^T,LD,TD refer to continuous portion of leading edge, side 

edge, trailing edge, discontinuity in leading edge, 
and discontinuity in trailing edge, respectively 

u,Z refer to upper and lower surfaces of airfoil, 

respectively 

Single or successive subscript coordinates indicate partial 
differentiation with respect to subscript variable. 

X ,y Indicates partial differentiation with respect to x 

and y, s and t times, respectively 


BASIC THEORY 

The analysis is based on the usual assumptions for thin airfoils 
in the linearized potential field. A solution for the disturbance 
parameters can thus be obtained by integrations of source and doublet 
singularities in the plane of the wing (z = O). The basic equation is 
(for example, references 9 and lO) 

fi(x,y,z) = - - 57^) + ^ (?) d’l (1) 

T 

where the function 0 is a solution of the linearized partial differ¬ 
ential equation for subsonic and supersonic flows 



In equation (l), r = V(x-e)2- [iiy-Ti)*^+z“J and the region t 
includes the entire z = 0 plane that can influence the point (x,y,z). 
At subsonic speeds, the factor K is equal to l/Zn. At supersonic 
speeds, the factor K is equal to l/m and only the finite part of 
the integral is used. It is important to note that the function can 
represent either the velocity potential or any of its derivatives, ai.d 
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D 

if all these values vanish sufficiently far ahead of the wing, the 

Q are solutions of equation (2). In case the^deriya- 
tive^f Ci becomes infinite at one or more points, the subsoitution 
of this derivative for in equation (2) depends on the condi o 
that the isolation of each singularity yields a finite integrand in 

the limit. 

SOURCE DISTRIBUTION FOR DOWH¥ASH FIELD 

A lifting wing of infinitesimal thickness is considered. The 
nerturbation velocity u vanishes everywhere in the z - 0 plane 
Except on the wing itself, and the perturbation velocity ^ 
everywhere in the z = 0 plane except on the wing and in the wake, 
ll e^ation (l), letting Q = w and noting that 


w - w, = 0 
u I 




° - 'Vi 


result in 


w(x,y,z) = - kJJI (B^u^-v^P de dq 


where T represents one surface of the wing and of the wake that can 
influence the point x,y,z. The quantity v^^ is a function of u 
through the irrotationality relation 


V = u. 


Therefore, 
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According to equation (3), the perturbation velocity w at any 
point in the flow field due to a lifting surface is determined by an 
integration of elementary source solutions with the distribution 
function given by (B^u^-v^). 

THE SEMI-INFINITE OBLIQUE WING WITH UNIFORM LOADING 

An important application of equation (3) can be obtained by con¬ 
sidering a semi-infinite oblique wing with uniform loading (fig. l). 
.The origin of the coordinate system is taken at the intersection of 
the leading and side edges, where both edges extend to infinity. For 
uniform loading, the term u^ in equation (3) vanishes everywhere over 
the region t except across the leading edge. In evaluating the dis¬ 
continuities in u or V that occur across the edges, a limiting 
procedure is used throughout the present analysis, which corresponds 
to the assumption of a linear variation in u or v across an 
infinitesimal strip of the edge. The distribution of v over the 
region t for this type of loading, as obtained from equation (4a), 
is shovm in figure 2, and v^ is seen to vanish everywhere except 
across the edges. In terms of equation (3), therefore, the downwash 
field for a semi-infinite oblique wing of uniform loading is obtsdned 
by means of a line integration of sources along the edges. 

It is subsequently shown that the downwash solutions for the semi- 
infinite line sources along the leading and side edges may be used by 
simple manipulation to obtain the downwash field for finite plan forms 
of uniform loading. 


The Semi-Infinite Oblique Leading Edge 

The semi-infinite oblique leading edge with origin at (0,0,0) con¬ 
tributes to u| and Vyj along the edge and, by virtue of the dis¬ 
continuity of the edge at the origin, contributes to v^ along the 
x-axis beginning at the origin and extending backwards to infinity. 

(See distribution of v in fig. 2.) The downwash contribution of 
the continuous portion of the leading edge is designated W]^; whereas 
the contribution of the edge discontinuity at the origin is designated 
^LD- 
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Contribution of continuous portion of leading edge . - The integra- 
■tion of equation (3) along the oblique leading edge is most conveniently 
performed in terms of an oblique system of coordinates, such that 


x' = x-iubV 
y’ = y-mx 


(5) 



It may be shown that the differential equation (2) is invariant under 
the change of variables x—►x', y-*«y', and z—*>z’ at both subsonic 
and supersonic speeds. If 0 (x,y,z) is a solution of equation (2), 
Q(x',y',z') is therefore also a solution. For other examples of the 
use of the oblique transformation in wing-theory problems, see refer¬ 
ences 11 to 13. 

In the oblique coordinate system, equation (3) evaluated along 
the wing leading edge becomes 


\^(x',y',z') = -K r ;^w^.ATi’d|' 

where 

r’ = B^(y'“+z'^) 


Noting that 


(Wz')t = 


B^u ^>-v^ t-mB^(u^^,-v^i) 




( 6 ) 


and along the leading edge for an infinitesimal strip of width Atj', 


-5.- Vj.-O 


u 


v^.An' = 


_o 

m 




T)^ * 0 
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results in 


;^(xSy’,z') = K 


Un VI-BV 


At subsonic speeds, the upper limit 4 2 - *•■'^*3, the integral is 

divergent. In actusil cases, however, this divergence does not present 
any difficulty because the construction of the finite leading edge by 
means of the superposition of two semi-infinite leading edges of 
opposite sign leads to the result that the infinite upper limit 
cancels. At supersonic speeds, i '2 is the position of the last 
source with Jfech aftercone including the point x',y’, 2 ’j that is. 


= x'-BVy‘^+z'2 


Integration of equation (7) yields the following expressions for 
the line source originating at (0,0): 


At subsonic speeds. 


where 


d 




B^ = 'sjl-V? 


In equation (8a), the term that etrises from substituting the upper 
limit of integration has been neglected because, as indicated pre¬ 
viously, it vanishes for leading edges of finite lengths. 

At supersonic speeds, for Bm < 1 (subsonic leading edge), 


■Wf (x,y,z) 


cosh 


-1 X' 
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for Bm > 1 (supersonic leading edge), equation (8b) becomes 


WL(x,y,z) = - ^ Vb ^®^-1 ^ fr 

d 


In equations (8) and in all subsequent ejcpressions, the positive 
value in a radical term must be preserved when extracting the root. 
For example, if 

y < 0 


then 


V? = V(-y)^ = -y = |y| 


For a leading edge normal to the flight direction, m = » and 
equations (8) yield the following expressions: 


At subsonic speeds^ 


^ ^I’^O ..-...-n ®iy 

WL(x,y,z) = — srnh ; 




At supersonic speeds , 


. X ^0 ^-1 -By 

WL(x,y,z) = - — cos -y== 


Vx2-b2z2 

It is shown by equations (8) that the downwash fields wj_^ con¬ 
tributed by the oblique and nornial leading edges are conical, or w 
is constant along radial lines emanating from the origin. 

Contribution of leading-edge discontinuity at origin . - Frc^in fig¬ 
ure 2, it can be seen that the leading edge cut-off at the origin 
results in 

V At) = - -2 (9i 

q m 


Along the semi-infinite side edge. \'Jhen equation (9a) is substituted 
into equation (3) (with u^ equal to zero along the side edge), there 

results 


2121 


WACA TN 2344 


11 





(9b) 


where ^2 position of the last source that can influence the 

point XfjfZ} that is^ at subsonic speeds, ^2 infinity, and at 

supersonic speeds, 

^2 = x-B^j^+z^ 


Integration of equation (9b) yields the following: 

At subsonic speeds, 

vj.p(x,y,z) = - ^ sinh-1 ^ ^ (lOa) 

where the term arising from the infinite upper limit of integration has 
been neglected because it vanishes for leading edges of finite length 
and 

ea = VIb^I (y2+22) 


At supersonic speeds. 


r » , _i X 

w (x,y,z) --cosh j— 

LD' ' jtm 5 a 


(10b) , 


It is shown by equations (lO) that the downwash contribution of the 
leading-edge discontinuity at the origin resialts entirely from the 
obliquity of the leading edge. Thus, if the leading edge is normal to 
the flight direction, this contribution vanishes. 


The Semi-Infinite Streamwise Side Edge 

Along the streamwise side edge, the quantity u«. vanishes; thus 
the source-distribution function along this edge is proportional only 
to the quantity v,^. 

The details of evaluating the contribution of the streamwise side 
edge to the downwash field are given in the appendix. The results for 
the semi-infinite left streamwise side edge are as follows: 
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At subsonic speeds. 


W3(x,y,z) 


f 1 , 

\x +*\jy?+i^ id y 


(11a) 


At supersonic speeds, 


ws(x,y,z) 


B^^py f 1 x\ 

« \x +Vx^-id^ 


(11b) 


Equations (ll) show that the downwash contribution of the streamwise 
side edge is independent of the obliquity of the leading edge. 




FINITE WINGS WITH UNIFORM LOADING 

The downwash field due to finite wings with uniform loadi^ can ^ 
be obtained by superimposing the fields due to a number of semi-infinite 
wings of the type considered in the previous section (fig. 1). This 
superposition is equivalent to superimposing the fields due to the 
source lines expressed by equations (8), (lO), and (11). 

A plan fom with curved edges requires an infinite number of 
source lines. If the edges are composed of straight-line segments, 
however, a finite number of source lines can be used to represent the 
plan form. 


Swept Wings with Streamwise Tips 

Let the downwash field due to a semi-infinite line source 
beginning at XQ,yQ be denoted by w^ let the subscripts 

L S, T, LD, and TD refer to the leading, side, and trailing edges 
and the leading- and trailing-edge discontinuities, respectively, all for 
the right half-wing; and let w refer to the effect of the source 
lines originating on the left half-wing. Then, the downwash field for 
a unifornly loaded wing, which is symmetrical with respect to the 
x-axis and which has uniformly swept leading edges and streamwise tips, 
is given by the following sum; 
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= 


(w+w) 


0,0 


(w-w) 


® -"L,LD 


(w+w) ■ (w+w) 

‘^rjO h 


+ct,*h 


T,TD 


(w+w) 



(v+v) 

h 


m 


+c 



(12a) 


In equation (l2a) and in the subsequent expressions, the upper and 
lower signs preceding a term refer to the right and left half-wings, 
respectively. 

The downwash fields due to the semi-infinite line sources indi¬ 
cated in equation (12a) can be obtained by simple manipulation of equa¬ 
tions (8), UO), and (11). The following transformations are made to 
obtain the effect of a semi-infinite line source originating at XQ,yQ: 


X in equations (8), (lO), and (ll) = x-Xq 
y in equations (8), (lO), and (ll) = i(y-yo) 


(12b) 


Equations (8) and (lO) can also be applied to a semi-infinite line 
source along the trailing edge by replacing m in these equations 
with the cotangent of the sweep angle of the trailing edge. 

At subsonic speeds, every point in the field is affected by all 
the terms in equation (12a); whereas at supersonic speeds, the point 
is affected only by those terms that refer to edges which lie within 
the Mach forecone from, the point under consideration. 


Rectangular Wings 

For the rectangular wing, the edge discontinuities LD and TD 
disappear and equation (l2a) becomes 
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WINGS WITH AEBITRARY LOADING 

If the wing loading is an arbitrary function of x and y, the 
integration in equation (3) for the downwash field is^ in general, 
required over the entire wing siirface and in the wake. An alternative 
•procedure in this case is the use of "correspondence formulas" as 
indicated in reference 8 by means of which the downwash field due to 
a variable loading may be expressed in terms of the downwash field due 
to a uniform loading with the addition of corrections for the edges of 
the plan form. 

In reference S, it is shown that if u on the wing is e^ressed 
in a series as a function of x and y in terms of the uniform^ 
prescribed velocity Uq and any term in the series is differentiated 
with respect to x and y, s and t times, respectively, such that 

(u) s t = b'iQ wing) (13a) 

X y 

where b is a consban't> ’then subject, to edge corrections^ there 
results 

('^n,k) s t = bWQ(x,y,z) (l3b) 

X y 

where w , and w^ refer to the downwash fields due to the term in 

n^ iC 0 

the series and to a uniform loading, respectively. The edges of finite 
plan forms may alter the given relation (l3a) so that edge corrections 
may be required for the relation e 3 q)ressed by equation (13b). At 
supersonic speeds, the relation given by equation (13b) thus applies 
to finite wings at all points in the flow field outside of the Mach 
aftercone from the edges that alter the given relation (13a) on the 
wing. Correspondence formulas for rectangular wings at subsonic and 
supersonic speeds are given in reference 8, table I. 
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ILLUSTRATIVE EXAMPLE FOR RECTANGULAR WING 
AT SUPERSONIC SPEEDS 

As an illustration of the method described herein, the downwash 
field is obtained at supersonic speeds for a rectangular wing with 
uniform loading and with linear chordwise variation in loading. The 
lettered regions in the subsequent discussion refer to figure 3. 


Uniform Loading 

For the rectangular wing with uniform loading at supersonic speed, 
equations (8e), (Ub), and (12c) are applicable. 

Region A . - In the region within the leading-edge I&ch cones and 
outside of the side-edge and trailing-edge Mach cones, 

(vQ)j^{x,y,z) = (wo,o+^0,o)l 



Region B . - In the region within the trailing-edge Mach cones and 
outside of the side-edge Mach cones, 

(Wo)B(x,y,z) = (Wq)^ - (w^^Q + 



(15) 


Region C . - Region C, which refers to points outside the trailing- 
edge Mach cones, is divided into two subregions depending on whether 
the point lies within one or both side-edge Mach cones. 

Region C 2 _. - If the point lies within one side-edge Mach cone. 
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(wo)ci(x>y>^) 


where 


and 


Region Cp 


(vo)c2(x^y^z) 


where 


and 


Region D » 
cones from the 
divided into f: 

Region 


= (wo)a - (wo,h)L - (^0,h)s 



(16) 


ea = BVya^+2^ 


ya = y-h 

- For points within both side-edge Mach cones > 


~ ^Vct " ^'^0,-hh ■ ^''o,-h^S 


BUq 


cos 


Byb 




+Byb 


x+ 


Vx^-lv,^ 


t 2 


(17) 


4 


yt = y+B 

- Region D refers to points that are always within Mach 
leading, side, and trailing edges. This region is sub- 
.ve regions, as shown in figure 3. 


- For points in D-j^, 


+ V 


c,0 c,0 T 


),r = (^n)n + 


O'C- 


0 


(18) 
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Region D 2 . - For points in Dg, 




■ <“o)c, " ^ 


COS 


-1 By, 


- - ^ 1 

^yix-o) 


By. 


a 


^x-c+yu-c)2-ea^ i 


x~c 
2 
a 


(19) 


Region Dj. - For points in D 


S-’ 




- (^o,.ak - (%-h)s 

~ . +By^ 

(_ Vx2-b2z2 


ih 


t20) 


Region D 4 . - For points in D 


4^ 




3 •*^1 

( 21 ) 


Region D 5 . - For points in D 


5^ 


= (^ok^ + (^c,.h)T + (Wc,.h) 


- (v ■> 

- (^ 0 )d 4 - 


COS 


-1 


By. 


yix-c) 


By, 


x-c 
2 


x-c+y(x.c)^-eb^ 


( 22 ) 
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Linear Chordwise Loading 


For linear chord-vdse loading, let 


u, = ^^0^ wing) 


(23a) 


TViP downmsh field for this type of 
where a and Uq are cons ^ • ^^^^.^ggpondence formulas of refer- 

loading can be obtained by obtained in the preceding 

ence 8 in conjunction ving at supersonic speed, 

section for the unifor y o -fields due to linear chordwise 


rr ^aSd T ^nffrt^ r~~ds"due to n^ear chordwise 

If w^ an Wq - _ t.hen from reference 8, 


If w^ ana w^ respectively, then from reference 8, 

«a=. cones, 


w-,(x,y,z) = a r WQ(^,y,z) d| 

•Jbz 


(23b) 


and for points within the trailing-edge Mach cones. 


W3^(x,y,z) ^ a 


wo(^^y>z) d^+ c(AWo)t(x 




(23c) 


.here (Awo)t 

0 ( 3 ,g 0 fox* iinifox*in. los^ing* 

« • A -Po-r the region within the leading-edge Mach cone and 

outsliSEviio-eage end tralllng-edge ifcch cones, eguntlons (23b) 
and (14) are applicable and there results 

(Wi)A(x,y^z) = -sBuq (x-Bz) 

T, • w Tn the region within the trailing-edge Mach cones and 

Region B* - 121 tne g __and (15) yield 

outside th4"side-edge Mach cones, equations (23c), (1 ), 

pc+Bz 

(Wi)B(x,y,s) - -sBno J " ° 

„ ■ r Tr, the region ulthln one or both side-edge ®ch cones 

and oiS&^ie'tSH4nS“Mocb cones, equation (23b) applies. 
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Il(x,ya,z)= X cos' 


By^ _i X _ -1 

^ - - -By. cosh""^ j- - Bz cos - / V" ? 

'^ 2 . 32^2 ^ UVx 


(26d) 


The solution for Ip is 


l2(x,ya,z) = | cosh‘l ^ 


r 2 _? 2 


(26e) 


Utilizing the solutions for and I2 given by equations (26d) 

and (26e), equation (26a) yields for points in region C within one side- 
edge Mach cone, 

(wi)c^(x,y,z) = - ^ [li(x,yg^,z) + Byg^l2(x,ya,z)] (26f) 

Region C 2 * - For points in region C within both side-edge Mach 
cones, equation (23b) may be written 




= a r W (i, 

Jnrr ^ 


= ("'l) C. + ~ 


y^z) 


'bz 


i L4'- 




(27a) 


where Vo(^jy>z) been obtained from equations (16) and (17). 

Utilizing the solutions for Ii and Ig given in equations (26d) 
and (26e), in which and are replaced by y^, and 1^, 

respectively, changes equation (27a) to yield 


(Wi)c2(^»y^2) = (Vc, 


rii(x,y^. 


z) + By, l2(x 




(27b) 
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Region D > - In the region vithin the trailing--edge Mach cone and 
within one or both side-edge Mach cones, equation (23c) applies* 

Region - For points in region the integral in equa¬ 
tion (23cy becomes 




+ aBu^ (x-c-Bz) 


(2£a) 


where WQ(C,y,z) has been obtained from equations (16) and (18). 

The quantity (Awq)^ in equation (23c) is evaluated for this case 
as the difference between (wq)d^ and (■Wo)c^i that is, from equa¬ 
tion (18), (AWq)^ = these results are combined, equa¬ 

tion (23c) yields 


(wi)£|^(x,y,z) = (Wj^)(,^ + aBu^ (x-Bz) 


(28b) 


Region D 2 . - For points in region D 2 , the integral in equa¬ 
tion (23c) becomes 


^Bz 


aBu 


y,z)de= (wp^^ + — 


r 

•-'c+Bz 


cos 


-1 


By, 




2„2 


d^ + 


B^z 


By. 


'c+e„ 


l-c 


ea' 


de 


(29a) 


where WQ(^,y,z) has been obtained from equations (16) and (19). The 

integrals in equation (29a) may be evaluated by utilizing the solutions 
for and I 2 as expressed by equations (26d) and (26e). In these 

solutions X is replaced by x-c. The quantity (AWq)^ in equa¬ 
tion (23c) is evaluated in this case as the difference between (w„) 

0 D 2 

and • When these results are combined, equation (23c) yields 
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(vi)D.Xx,y,z) = (w^) + Il(x-c,yg^,2) 4 - By^^IgCx-c^yg^,z) + 


-1 _i-D, / 1 x~c 

c cos — _ - ^ +Bya ( . . 

‘V(x-c)^-B^z^ \ X-C+ > 


Region D3. - For points in region D3, equation (23c) yields 


(29b) 


(w.)nJx,y,z) = (w^)^ + 




Hr. 


i_ - -i.Ue 

2 S 2 t. 21 




(30a) 


where WQ(E,y,z) has been obtained from equations (18) and (20) and the 

term (Aw-)^, is included in (w ) . The integrals in equation (30a) 

U 1 -L 

may be evaluated by utilizing the solutions for and Ig as 

expressed by equations (26d) and (26e). In these solutions, y^ is 
replaced by y^; thus 

aBuQ r “I 

(wi)D^(x,y,z) = (^3^)0^ + -7- [_Ii(x,y^j,z) + By^ l2(x,yb,z^ (30b) 

Region D4. - For points in region D4, equation (23c) yields 

(wi)j) (x,y,z) = (wp + (vp - (w^) 


where Wf>(C,y,z) has been obtained from equations (19) and (21). 


(31) 
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Region D 5 . - For points in region equation (23c) yields 


(wj)D_(x,y,z) = (w ) 


aBu, 


0 ^ 


1'D4 n 


COS 


-1 


'c+Bz 




dC + 


Byi 


i-c 




di + 


cos 


-1 




y(x-c) 


+Byb 


x-c 

2 


X-C+ 




(32a) 


where WQ(i,y,z) has been obtained from equation (22). The integrals 
in equation (32a) may be evaluated by utilizing the solutions for and 
I 2 as expressed by equations (26d) and (26e). In these solutions, 

X is replaced by x-c, and y^^ and are replaced by and 

respectively. Thus 


(Wi)j2^(x,y,z) = (wt) 


aBu 


0 




1 


cos 


-1 


Il(x-c,y^,z) + By^, l2(x-c,yyj,z) + 

Byb 


V(x-c)2-B^z^ 


+By> 


X-C 


x-c+y(x-c)^-|^^ Jj 

(32b) 


Levis Flight Propulsion Laboratory^ 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, January 22, 1951. 
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APPENDIX 

DERIVATION OF DOWNWASH FIELD DUE TO SEMI-INFINITE 
LEFT STREAMVniSE SIDE EDGE 

The subsequent derivation for the dovnvash field refers to the 
left streamvfise side edge of the uniformly loaded semi-infinite oblique 
wing, as shown in figure 1. 

The downwash field contributed by the side edge may be obtained by 
means of equation (3), in which the integral is evaluated along the side 
edge. As noted in the text, the quantity u^ is zero along the stream- 

wise side edge (of width An); therefore. 

Across "the side edge ihere is an. abrupi increase in u* This 
increase is assumed iso occur over an infinitesimal width € across the 
side edge. Then according to equation (4a), 

v(x,0,0) = r u;r|‘ii= ~~ ^ 

Jo ^ 

y =s £ ^ however, u.^^ due to the side edge vanishes and therefore 

v(x,€,0) = 0 (•^^) 

The distribution of v in the vicinity of the side edge (y = O) 
is shown in the following sketch; 



y- 
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It is assumed that the abrupt changes in v indicated in the 
sketch occur over the infinitesimal distances At), so that along the 
line (CjO,0), 




Substituting these vedues for into equation (Al) results in 


Wg(x,y,z) 





where the two upper limits refer to subsonic and supersonic speeds, 
respectively. Integration of equation (A6) yields 
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vrnen the limits are substituted into equation (A7) and e is 
made to approach zero, the following results are obtained: 

At subsonic speeds. 


•ws(x,y,z) 


2n 



(A8) 


where 


= V|b^| (y^+z^) 


At supersonic speeds, 


V3(x,y,z) 



(A9) 
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A method utilizing source singularities is pre- A method utilizing source singularities is pre¬ 
sented for obtaining the linearized dovnwash field due sented for obtaining the linearized downwash field due 

to lifting wings at subsonic and supersonic speeds. The to lifting wings at subsonic and supersonic speeds. The 

method is applied to derive generalized formulas for the method is applied to derive generalized formulas for the 

downwash field due to uniformly loaded swept and rectan- downwash field due to uniformly loaded swept and rectan- 
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gular wings at subsonic and supersonic speeds. The gular wings at subsonic and supersonic speeds. The 

utilization of these formulas to obtain the downwash due utilization of these formulas to obtain the downwash due 

to wings of arbitrary loading is indicated. to wings of arbitrary loading is indicated. 







